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Abstract 

We compare the Infrared Dirac-Born-Infeld (IR DEI) brane inflation model to observations 
using a Bayesian analysis. The current data cannot distinguish it from the ACDM model, 
but is able to give interesting constraints on various microscopic parameters including the 
mass of the brane moduli potential, the fundamental string scale, the charge or warp factor 
of throats, and the number of the mobile branes. We quantify some distinctive testable 
predictions with stringy signatures, such as the large non-Gaussianity, and the large, but 
regional, running of the spectral index. These results illustrate how we may be able to probe 
aspects of string theory using cosmological observations. 
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1 Introduction 



1.1 Experiments and theories 

An ongoing and forthcoming array of experiments (e.g. WMAP, [l], SDSS [2], [3], HI [5] , SNLS 
[6], ACBAR [7], Planck [H], ACT [9], Spider [lO]) is measuring the cosmic microwave back- 
ground (CMB) and the large scale structure of the universe with unprecedented precision. 
This provides exciting opportunities to reveal the nature of the early universe and the un- 
derlying fundamental theories. The leading theoretical candidate for creating the initial 
conditions of our universe is inflation [HI [121 US]- However inflation remains a paradigm, 
which can be implemented by a variety of models underpinned by differing microphysical 
constructions; as the constraints from data tighten, there is the hope that we might identify 
the specific scenario that describes our universe. With the natural ingredients of such model- 
building being supergravity and string theory, the process of better measuring the properties 
of the early universe is also a process of understanding better the theory of quantum gravity. 

In contrast to recent debates on the predictivity of the string theory landscape, here 
we use a more conventional approach to investigate the predictivity of string theory by 
studying the properties and exploring the dynamics of our own vacuum. We first scan the 
parameter space of inflationary models subject only to the requirement that they provide 
enough inflationary e-folds to solve the flatness and horizon problems. This is because the 
natural creation of a homogeneous and isotropic universe is the leading problem that we want 
to solve, and is perhaps the most attractive feature of the inflationary paradigm. After that, 
we study the observational consequences of all the viable parameter spaces, with the goal 
of looking for distinctive signatures. Some of these can be compared with observations and 
used to narrow down the parameter space. Despite of the vastness of all possible vacua in 
the string landscape, this process can be rather effective since certain observational features 
rely on distinctive dynamics. As we will see, such dynamics can either be field-theoretic with 
strong motivations from string theory, or completely stringy in nature. 

There are many candidate observable signatures in inflationary models. The most generic 
ones are the amplitude of the primordial power spectrum and its spectral index. Since most 
viable models built from a fundamental theory have adjustable parameters to fit these two 
observables, this leaves a large number of viable models that are consistent with the data, 
and even leaves the nature of the inflaton field ambiguous. In principle. Nature is not 
obligated to provide more information within our experimental abilities, and indeed there is 
no evidence for further parameters required to describe the current data. But anticipating 
her generosity, possible distinctive observables that might be measurable in the future include 
the scale-dependence ( "running" ) of the spectral index, departures from Gaussianity of the 
primordial fluctuations, a tensor contribution to the primordial power spectrum, and cosmic 
strings. These will be crucial to successfully carry out the program that we have outlined. 

With the rapidly improving quality of cosmological data, it will become increasingly in- 
teresting to implement the above program by comparing specific models to data, starting 
directly from microscopic parameters of theories. Modern cosmological data analyses make 
use of the powerful method known as Markov Chain Monte Carlo (MCMC) to implement 
the comparison to data, providing an efficient way of estimating posterior distributions of 
the microscopic parameters. However, in practice, when directly using microscopic param- 
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eters as MCMC parameters, highly non-hnear relationships between the parameters and 
observables may introduce severe obstacles for MCMC to efficiently search the parameter 
space. Therefore, a reparameterization according to the specific nature of the model often 
becomes necessary. So instead of a straightforward exercise, implementing MCMC becomes 
a rather interesting model-dependent art. It is also a purpose of this paper to use an ex- 
ample to illustrate this process and extract certain model-independent procedures of such 
reparameterization which may be of more general interest. 

1.2 Brane inflation 

The inflationary models that we study in this paper belong to the brane inflation scenario 
proposed by Dvali and Tye [HI [15] Q We are interested in these models precisely because 
they can give rise to a large number of distinctive observational signatures. This happens 
even in the simplest scenarios that provide inflation. One of the most important reasons 
that makes it possible is that brane inflation can be achieved via two different mechanisms, 
namely slow-roll and Dirac-Born-Infeld (DBI) inflation. 

The original models of brane inflation [TU [20], [211 [22] are slow-roll inflationary models 
[T2l [13] . where branes and anti-branes slowly approach each other in a flat potential. A 
model that uses this mechanism in the framework of the string theory flux compactification 
[2B] is studied by Kachru, Kallosh, Linde, Maldacena, McAllister and Trivedi (KKLMMT) 
|22j . As in the F-term inflation models in supergravity [23], it is found that the generic 
shape of the potential is too steep to achieve the slow-roll inflation, in this case due to 
the moduli stabilization. Again, similar to those supergravity models, it is possible that 
several contributions to the potential manage to cancel to a certain precision so that the 
potential becomes sufficiently flat. There are effective parameters in the model controlling 
the inflaton mass that can be adjusted to fit the observed spectral index [25]. The running 
of the spectral index, non-Gaussianities and tensor modes are all too small to be observed 
in the near future]^ 

Another inflationary mechanism that is so far uniquely found in brane inflation is the 
DBI inflation [521 [331 [331 [35]- ^BI inflation, the rolling velocity of inflaton branes is 
not determined by the shape of the potential but by the speed-limit of the warped internal 
space. Such warped spaces are naturally present in the extra dimensions due to fluxes used 
to stabilize the string compactification |36j . 

The ffist model that uses such a mechanism is that of Silverstein, Tong and Alishahiha 
(STA) [521 [33]. In this model, as the branes roll into a throat from the UV side of the 
warped space under a quadratic potential, its velocity gets restricted by the large warping 
in the IR side of the warped space. However, instead of having a potential with a generic 
mass term, a rather steep potential, characterized by a large inflaton mass, is required to 
achieve this UV DBI inflation. The reason is that, when the branes enter from the UV side 
of the warped space in the GKP-type warped compactification [3S], the energy provided 

^For recent reviews on other types of string inflation models, see Ref. [HI [171 [TH [19] . 

^Some observables become measurable if there are sharp features in the potential [26l [27l [281 [29] [30| . In 
addition, there are other important observational possibilities of brane inflation - cosmic strings and those 
related to reheating (ST] [15] - which apply to both slow-roll and DBI inflation. 
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by the antibranes sitting at the IR side is not large enough to drive DBI inflation even if 
there is the speed-limit, since the antibrane tension has been warped down correspondingly. 
Therefore an extra, steep, potential has to be added to raise the inflationary energy. In 
addition, embedded in the same warped compactification, the model generates large non- 
Gaussianities that exceed the experimental bound [251 [37], as well as excessive probe brane 
backreaction which we will address in Sec. 12. 1[ This is because in this model, the levels 
of non-Gaussianity and probe brane backreaction sensitively depend on the inflaton value, 
and it is viable only if the inflaton field is of (super-)Planckian size. However, the range of 
the inflaton field is restricted by some geometric conditions of the compactification and is 
sub-Planckian [381 [39l [iQl [25l [371 . 

To fully make use of the speed-limit of the warped space, it is better to make the branes 
roll out from the IR end, and use antibranes in other throats to provide the inflationary en- 
ergy. In this way the speed-limit of the branes and the inflationary energy become relatively 
independent of each other, leaving a rather flexible shape of the inflaton potential which has 
been the main problem of model-building. This is the model proposed in Ref. [Ml [35] . It 
can be generically realized in the multi-throat brane inflation scenario [34] . 

It happens that in this IR DBI inflation model, the large non-Gaussianities can also be 
small enough to satisfy the current observational bound [38] . This is partly because no matter 
how small the warp-factor (and consequently, how big the non-Gaussianity) the branes begin 
with, the level of non-Gaussianity decreases as the branes roll out and approaches its minimal 
value at the end of the inflation. Therefore, in the segment of the warped space traversed 
during the last 60 e-folds, the level of non-Gaussianity is among the smallest in the entire 
DBI inflation trajectory. Moreover, the geometric conditions that put a tight constraint on 
the STA model are automatically satisfied in the IR DBI model and has no effect on the 
non-Gaussianities. 

Besides providing a speed-limit to the inflaton, another important property of warped 
space is the reduction of the local fundamental string scale [H]. This turns out to have 
important consequences on density perturbations in DBI inflationary models. During the 
epoch when the string scale is red-shifted below the Hubble parameter, the quantum fluc- 
tuations on the inflaton branes become stringyH The density perturbations are no longer 
fully described by the usual field theory approximation, and acquire distinctive stringy sig- 
natures. In the IR DBI model, this stringy phase corresponds to earlier inflationary e-folds, 
and therefore larger scales in the sky. It is estimated that such a phase transition will give 
rise to a large transient (regional) running of the spectral index [551 ISSj- In this paper, we 
make this prediction more quantitative and compare it to observations. 

1.3 Outline 

Following the strategy that we outlined in Sec. 11.11 in this paper, we first summarize the 
overall features of brane inflation using phase diagrams that describe the parameter spaces 

■^Notice that such a stringy phase only happens in the inflaton sector, which is the deep IR side of a warped 
space with energy density of order H'^, so it does not backrcact significantly on the Hubble expansion. We 
also note that such a stringy phase will backrcact on the IR side of the warped geometry, but it is estimated 
that this still leaves a large enough portion of the geometry for DBI infiation to take place [35l [HI] ■ We will 
discuss this more in Sec. l2.2l fc l3.2l 
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spanned by both inflationary mechanisms, i.e. slow-roll vs. DBI (Sec. [2]), reviewing the key 
observational predictions in the different parts of the parameter space. 

The main focus of this paper is to compare the IR DBI brane inflation model to ob- 
servations (Sec. [3]). We derive analytical and numerical model predictions for the shape of 
the power spectrum, non-Gaussianity, and tensor modes, giving a quantitative estimate of 
the effect of the Hubble-expansion-induced stringy phase transition on density perturbations 
(Appendix El). 

We then proceed to compare these results to the observational data from cosmic mi- 
crowave background and large scale structure (Sec. Hj). We outline how such a comparison 
should be generally implemented using MCMC. The current data give a number of inter- 
esting constraints on the microscopic parameters of the model (Sec. [5]), including the mass 
of the brane moduli potential, the fundamental string scale, the charge or warp factor of 
throats, and the number of the mobile branes. We also quantify some distinctive observable 
signatures of this model, such as the level of the non-Gaussianity and the running of the spec- 
tral index. We discuss how the latter is observationally different from two other cases that 
may also give large running spectral index: slow-roll inflation with mild features on potential 
(Appendix [B]), and slow-roll or DBI inflation with a non-Bunch-Davies vacuum (Appendix 
ICl) . These results illustrate how string theory can make testable predictions which might be 
subject to observational constraints. 

For convenience, all the variables used in this paper are summarized in Table [TJ 

2 Phase diagrams of brane inflation 

As mentioned in the introduction, a useful approach to study inflationary models is to first 
scan through as large a parameter space as possible with the requirement of a sufficient 
number of inflationary e-folds. Then we can work out the observable predictions (such as 
density perturbations) in different regimes, and compare them with the data to narrow down 
the parameter space. 

Bearing this in mind, in this section we will study the parameter space in brane inflation 
models that can provide enough inflationary e-folds. In this paper we choose the repre- 
sentative examples in which D3-branes move along the radial direction of a throat with an 
approximate AdS geometry in type IIB flux compactification. 

In the case of the flat 4-d space-time and non-compact extra dimensions, D3-branes move 
freely in the throats. However realistic inflation models in warped compactification requires 
an inflationary 4-d space-time with a Hubble parameter H and stabilized compact extra 
dimensions. In this case the moduli space of branes are lifted and receive potentials with 
masses of order H. This is only the generic expectation - the details of the potential profiles 
are environmental, depending on various ingredients (such as fluxes and other branes) present 
in speciflc string compactiflcation models. 

In this paper we choose to simply parametrize such unknown mass terms in the potentials, 
hoping it can provide a bridge between the bottom-up observational data-fltting and top- 
down string theory calculations. 
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Variable Description 



Notes 



Mpi 


4d reduced Planck mass 


Mpi = (8^G)-i/2 ^ 2.4 X lO^^GeV 


nis 


Mass scale of fundamental strings 




9s 


String coupling 


5s < 1 


T3 


D3-brane tension 


Eq. D 


R 


Length scale of warped throat 


Eq. 


M, K 


Flux numbers in warped throat 


Integers 


UA 


Number of antibranes in A-throats 




Ub 


Number of branes (inflatons) in B-throat 




Nb 


Effective charge of B-throat 


Nb = qbMK 


aB 


Multiplicative factor from orbifolding 


~ 1 in data analysis 


\b 


As = nBNB/2'K^ 


Eq. ^ 


r 


Radial coordinate of throats 




4> 


Canonical inflaton field 


(j) = r^riT^ 


hA 


Minimum warp factor of A-throat 




HbU) 


Warp factor at location (j) in B-throat 


h = r/R = SR/^fX 


P 


Characterization of shape of potential 


Inflaton mass vn? = (3H'^ 


7 


Lorentz factor of inflaton 






Sound speed in 4d 


c, = 1/7 for DBI inflation 


Vo 


Inflationary energy density 




Ne 


Number of e-folds to the end of inflation 






Number of e-folds to the end of IR DBI inflation 




Nfl 


Total e-folds of non-relativistic roll inflation 


Typically fast-roll 


kc 


Critical scale of the stringy phase transition 


Eq. (j4.38|) 


Nc 


Critical DBI e-fold at kc 


Eq. (jA.7[) 


P{k) 


Power spectrum 






Tensor to scalar ratio 




feq 
J NL 


Estimator of the non-Gaussianity 


Equilateral shape 



Table 1: Description of variables. In the text, subscripts A and B are frequently added to some 
of the variables, referring to the quantities of the A- or B-throat. 



2.1 UV models 

In this and the next subsection, we draw the phase diagrams of brane inflation in terms of 
two parameters: the inflaton position and the mass of the inflaton moduli potential m. 
We use these diagrams to show the conditions under which different inflationary mechanisms 
happen. 

We first consider the UV models, the phase diagram for which is shown in Fig. [H In the 
UV models, the branes are started from the UV side of a throat (denoted as the A-throat) 
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and attracted to the IR end by the moduh potential or Coulomb potential from antibranes, 



V{<P) = H + ^mV' + V^Coulomb(0) , (2.1) 

where Vq = 2nAh\T^ are provided by ua antibranes at the end of the throat, which eventually 
get annihilated by the same number of branes. The warped geometry is 

ds^ oc -^dsl + , (2.2) 

where the ds^ is the 4-d space-time metric, and 

Xa = tiaTsR^ = uaNaI {2ti^) (2.3) 

where Na and Ra are the effective charge and characteristic length scale of the warped 
space, respectively, ua is the number of inflaton branes. Note that Na may include the 
multiplication factor a a from orbifolding on the original D3-charge Nqai Na = aAN^A- The 
following relations are also useful: 

4 

R\ = ^T^gsNAm-' , T3 = , (2.4) 

where = is the string mass scale. (Later we will also use the same definitions for 

other throats with corresponding subscripts A 01 B.) 

If the flatness of the potential V{(j)) satisfies the slow-roll conditions, the branes can 
slowly roll non-relativistically and the kinetic term of the brane DBI action reduces to the 
minimal non-relativistic form. To indicate this condition in the phase diagram, we draw a 
curve of riy = MpyV"/V = 1 which is 

0^ = 2M^, - ^ , (2.5) 



m 



where we neglected the Coulomb potential term for simplicity. Eq. (12. 5p corresponds to the 
solid line stretching from the lower- left (at m = ^/Vo/Mpl) to upper- right (at (p = \/2Mpi) 
in Fig. [H Inclusion of the Coulomb term will cause a slight deformation at the lower-left 
corner of this curve and will not affect our conclusion. The shaded region above and to 
the left of this curve has rjy < 1, and corresponds to the slow-roll inflation phase; here the 
condition < 1 is always weaker. The lower-left shaded region is the small-field slow-roll 
region where the potential is dominated by the constant Vq. As increases toward the upper 
part of the shaded region, it corresponds to large-field slow-roll inflation where the potential 
is dominated by m^0^ term. Note that in this region is trans-Planckian, > ^/2Mpl. 

Outside of this slow-roll region, naively the inflaton will roll down the potential very fast 
and make inflation impossible. However because of the presence of the warped space, the 
velocity of the inflaton is bounded by the warped speed of light and therefore cannot be 
arbitrarily increased. Silverstein and Tong [32] show that, if 



m > Mpl/^/XA , (2.6) 
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v^/Mpi 



> m 



Figure 1: The inflation phase diagram for UV models. The shaded regions correspond to parameter 
space that can give rise to inflation. The darker the region is, the larger e-folds it can provide. "S.R." 
stands for slow-roll inflation; "DBI" stands for DBI inflation. The arrows indicate the starting point 
and rolling direction of the inflaton. In brane inflation, inflatons have to stay below the horizontal 
solid line (at 0^ = Ra^uaT^)'-, the two vertical lines (at m = \/Vq/M-p\ and m = M-px/^/Xa) 
are widely separated. The curve stretching from m = ^/Vq/M-pi to (pA = V^M-p\ corresponds to 
r/y = /3/3 = 1. See text for discussion. 



the inflaton enters another phase of inflation, namely DBI inflation, in which the full form 
the DBI kinetic term has to be taken into account. We indicate this condition by the solid 
vertical line (at m = Mpi/V^I) in Fig. [T]and the DBI inflation phase by the shaded region 
to the right. 

There are two possible regions where these two phases merge onto each other: when the 
inflaton in the DBI phase starts from a Planckian value (the upper-right corner of Fig. [T]), 
the inflation can go continuously from slow-roll to DBI; when a/T^/^pi ~ ^pi/a/Aa so the 
two vertical lines in Fig. [1] become very close to each other or even switch places, the inflaton 
around this border will trigger an inflationary phase that lies in-between the slow-roll and 
DBI regimes. These are the "intermediate regions" studied in Refs. |^ [25]. 

However, so far we have been discussing an effective field theory description where one is 
allowed to independently choose the throat charge Na, fundamental string mass and the 
inflaton field range 0. In a realistic string compactification like the generic multi-throat flux 
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compactification in type IIB string theory [SB], various throats are glued to a bulk. There 
are several fairly model-independent geometric conditions that should be imposed to be con- 
sistent with the brane inflation setup that we have in mind. In this type of compactification, 
the Planck mass is obtained by an integration throughout the compact space, 

Ml, ~ g:'m% , (2.7) 

where Vq is the total volume of the compactification (after modding out the possible orbifold 
effect) and its dominant contribution comes from the bulk and the UV regions of throats. 
The throats are glued to the bulk and their sizes R are restricted, 

— <Ve. (2.8) 

aA 

The volume of the throat is divided by a a in case of orbifolding. The inflaton brane separation 
is restricted by 



0<i?AV^^, or 0<Lv/^, (2.9) 

for Ha branes moving in the throat or the bulk respectively. Here L is the size of the 
bulk in a certain dimension. These conditions have been used in various contexts in brane 
inflation, to constrain slow-roll models in the bulk [20l [22], random walk eternal inflation 
[39] . the UV DBI model through the non-Gaussianities [381 SOI [251 [37], ^-^d the tensor mode 
For our phase diagram these imply 

^/""^^ <^«1, (2.10) 



Mpi / ~ NA/aA 

and 



RaV^_ < vg^ 1 _ (2.11) 



PI 



^/ NA/aA 



Thus the two vertical solid lines (at m = a/Vq/Mpi and m = Mpi/ a/Aa respectively) in Fig.[T] 
should be widely separated; the inflaton can only move below the horizontal solid line (at 
= Ra yJuAT'i) which is well below = -\/2^pi- This excludes the large field models where 
> Mpi, and opens up a wide region in the middle of the parameter space where there is 
no inflation. This is also why random-walk eternal inflation in KKLMMT model within the 
throat is excluded in [32], the tensor mode in brane inflation is unobservable [ID], and the 



"intermediate UV models" are ruled out in [25] o The earlier statement that the antibrane 
tension alone is not large enough to drive DBI inflation in UV models can be justified as 
well by estimating i^anti^^ ~ Vra^Ta^R^/i^/Mpi -C 1, where i^anti is the contribution from 



^ Note that in the bulk, the geometric conditions described here alone arc not enough to restrict the 
scalar field to be sub-Planckian. For example, for a toroidal compactification with Li = UmJ^ (i = 1, . . . , 6), 
consider an irregular shape li > g^^h • • ■ ^5- Additional consistency requirements are necessary, e.g. how to 
maintain the shape of the potential over Acj) of Planckian size while it is expected to vary over the string 
scale. 
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antibrane tension and At ^ RAh']^ is the time scale that the branes spend travehng down 
the throat. 

Having considered this inflationary phase diagram, we can now restrict the parameter 
space by comparing it to observations. In the UV model, we saw from the above discussions 
that there is a clean separation of slow-roll and DBI inflationary phases after the geometric 
conditions f l2.7l) -( l2l9l) are applied. The slow-roll region is the KKLMMT model |22] and is 
compatible with the current observations [251 H6| HTj. In this model, the infiaton mass may 
be adjusted to fit the spectral index. The running of spectral index and non-Gaussianities 
are unobservable if the potential is featureless. The tensor mode is also unobservable. The 
DBI region is the STA model [321 [33]. It predicts large non-Gaussianities with the estimator 

I/- I - 1.3^ , (2.12) 

(where p = m/ {\/QM-p\/ \/\a) ^1,) together with a possibly-observable tensor mode r ?s 
5/(p/fiVL) [33]. From the constraints ^B) and fl2TTD . one can see that, for one brane ua = 
IJnlI ^ V^i^A/ o-aY which cannot fit the observations [25], [381 HO]- One way to increase 
the field range is to increase the number of the infiaton branes. Here we emphasize another 
constraint coming from the relativistic probe brane backreactions discussed in [221 [S2]- In 
order to treat the mobile infiaton branes as probes of the warped background, Na/qa 3> uaI 
is required. Namely, the energy scale of the mobile branes cannot exceed the source of 
the warped background. On the other hand, combining (12. 9p . (12. lip and (12.120 . we have 
I/at^I > 1.3p^A^^/(a^?7,^). Using the relation I/atl] ~ 0.327^, these two requirements lead 
to p <^ 0.5, which is a contradiction. Note that this conclusion is independent of the value 
of ua and before any comparison with data is madeEl However we should note that this 
inconsistency appears when the STA model is embedded in the warped compactification 
of the GKP type, so it remains a viable field-theoretic model, and looking for other UV 
embeddings becomes an interesting question. 



2.2 IR models 

In the IR models, branes are started from the IR side of a throat (denoted as the B-throat) 
and roll toward the UV side under the moduli potential 

^(0) = K)-^mV- (2.13) 

The origin of cf) is at the tip of the throat (12. 2p . which can be realized for example if the tip 
of the throat is an orbifold fixed point. The Coulomb attraction from antibranes in other 
throats is neglected here unless / H'^ is very small. 

The IR model can arise in the following scenario [M] (illustrated in Fig. [2]): At the be- 
ginning, antibranes are naturally attracted to and settle down at the end of various throats 
induced by fluxes. However they are semi-stable at most, and will eventually annihilate 

■'Considering wrapped branes [48[ 149] effectively interprets ua in a different way, so it should be subject 
to the same conclusion discussed here. 
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• Antibrane o Brane 

S S 




Figure 2: Multi-throat brane inflation scenario. In the first figure, antibranes are settled down in 
throats. In the second figure, in some throats antibranes annihilate fluxes and generate branes. For 
a throat with tachyonic brane moduli, branes fall out and settle down somewhere else, triggering 
either IR or UV models of brane inflation. 

against some fluxes [50]. The end products are many branes. As mentioned, unlike an- 
tibranes, branes experience no potential if the extra dimensions are not compactly stabilized. 
But for realistic inflation models in string compactification, their moduli space is lifted. If 
the mass term is tachyonic as in (12.131) for a B-throat, these liberated branes will roll out. 
The inflationary energy Vq is provided by longer-living antibranes in other throats (denoted 
as the A-throats) or in the bulk. Shorter A-throats give more dominant contributions to Vq. 
These antibranes eventually get annihilated by some of the infiaton branes. The annihilation 
products will naturally heat low mass-scale sectors in case of tunnelling reheating [51], such 
as branes residing in very long throats or in a large bulk. 

In the absence of the warped space, the |?7v| = M^J = 1 line is 

02 = -2M2i + ^, (2.14) 

which corresponds to the vertical line (at m = ^/Vo/Mpl) in Fig. [3l Slow-roll inflation occurs 
when |?7y| < 1 which is the region to the left of this line. 

In the presence of the warped space, DBI inflation can be triggered even if the slow-roll 
condition is not satisfled. Ref. [SU [35] show that, for jS = / H"^ ~ \riy\ > 1, DBI inflation 
happens if 

< HRl^/^s . (2.15) 

Here, an important difference from the STA model condition (12.61) is that, because the 
inflationary energy Vq is provided by antibranes in other throats instead of the moduli 
potential itself, the shape of the potential that can achieve inflation becomes rather flexible. 
Of the special interest is the generic case ~ H^. Even for (3 < 1 when slow-roll inflation 
is possible, the speed-limit provided by the warped space cannot be neglected if the inflatons 
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Figure 3: The inflation phase diagram for IR models. The notation used here is the same as in 
Fig. [TJ The vertical line at m = y/\^/Mpi corresponds to \7]v\ = /3/3 = 1. The unshaded region 
may support a certain amount of non-relativistic fast-roll inflation. 



are started from the region [35] 



< 



. (2.16) 



This implies that the DBI and slow-roll phases can be smoothly connected by an intermediate 
region in this corner of the parameter space. Overall, the DBI inflation phase stays below 
the horizontal curve stretching from the origin to = HR^^/ubT^ in Fig. [3l The DBI 
region always stays well below the maximum inflaton extension (the horizontal solid line at 



RbV'^'bTs), since 



ffi?B<^==«l, (2.17) 



RbV^bTs ^ ^/NbJq^ 



where 



Vo = 2nAh\T^, (2.18) 



and (12.71) (12.81) are used, ua being the number of antibranes that get annihilated. We see 
that IR DBI inflation is completed in a very small region at the tip of the B-throat. Unlike 
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the UV DBI phase, the condition (12.91) is automatically satisfied. This also justifies the small 
field expansion in fl2.13p . 

We have treated the Hubble parameter H as a. constant. This can be verified using the 
geometric constraints (12. 7p and (12. 8p . because in the IR model the potential drop during 
infiation, estimated very conservatively, satisfies 



(2.19) 



Vo ~ 2Vo NB/aB 
As long as 

(3 < NB/iaBUB) , (2.20) 

the infiationary energy is approximately a constant. 

Therefore we see that in the IR models, infiation can occur for a large range of the mass 
parameter, 

< m^/H^ < Nb , (2.21) 

around the generically expected magnitude m ^ H. The requirement is to start the infiatons 
from a small enough 0^. In terms of the flux compactification this is easy to achieve, because 
the minimum warp factor is given by the flux numbers M and K in an exponential form [36] 

/i^in ~ exp(-27rJ^/3M(7,) . (2.22) 

We emphasize that non-trivial constraints come from various back-reactions that cut off the 
IR regions of a throat. These include the back-reaction from the 4-d inflationary background 
[35l[5T] and the back-reaction from the relativistic inflaton branes [321 [52]. The former is gen- 
erally (for /3 ~ 1) more important than or comparable to the latter depending on the number 
of inflaton branes. It is estimated [521 EI] to cut off the throat at ~ HR^^ubT^/ ^/T^i 

This determines the maximum number of e-folds achievable by the DBI inflationary phase 
rDBi 



-^tot ~ V^B- A more detailed understanding of this backreaction is important. 



'tot 

It is worth pointing out that, in terms of model building, the most important difference 
between the IR and UV DBI models is not whether branes are started from the IR or 
UV side of a warped space. It is the independence between the inflaton speed-limit and 
the inflationary energy, which allows a flexible shape of potential. We have seen that this 
naturally happens when branes are moving out of a throat, with inflationary energy provided 
by antibranes in other throats. Just in terms of fleld theory, even in the UV models, such 
an independence can be achieved by demanding the constant term Vq in the potential (12. ip 
to be independent of the A-throat warp factor, for example by a hybrid of a different fleld 
around = to suddenly end inflation. The question is then how to realize it naturally in 
string models. 

Having considered the phase diagram, we would like to flrst restrict the parameter space 
by comparing the predictions of the model with observational data, and then make predic- 
tions for future observations. This will be the main focus for the rest of the paper. We close 
this subsection with a few comments on the setup of the model. 
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^This constraint is equally important for UV DBI models. 
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Firstly, as we demonstrated in Fig. [H for UV models there is no inflation around rjv ^ 1 
and a large parameter space beyond that. So for IR models, it is reasonable to consider the 
simplest case where, after branes come out of the B-throat, there is no signiflcant amount 
of additional inflationary e-folds if they roll through the bulk or enter another A-throat 
to annihilate antibranes there. This simplest possibility represents a fairly generic class of 
models. 

Secondly, more realistic throats such as the Klebanov-Strassler throat [53] have a scale- 
dependent charge. The characteristic scale R decreases slowly towards the tip of the throat. 
Especially the geometry around the tip region will be signiflcantly different from (12. 2p . For 
UV models such modiflcations can be important because the tip of the throat is the region 
around which the last 60 e-folds of inflation happens [54l[55]. For IR models, the situation is 
opposite. The last 60 e-folds of inflation happens away from the bottom of a throat because 
generally the total e-folds is more than 60. Furthermore the relevant fleld range is very 
small. Therefore, under these conditions (which will be made more precise in Sec. 13.11) . we 
can ignore both the deformation of the throat geometry and the running of the throat charge, 
and approximate the metric as (12. 2p with the constant A (or R) being the effective value at 
the relevant (or r). 

Thirdly, as we have mentioned, the realistic IR case almost always involves multiple 
inflaton branes. It is interesting to see whether the non-Abelian action plays an important 
role [56]. In our scenario, after the mobile branes are created in the IR end of the throat, they 
all have approximately the same radial coordinates and roll in the radial direction. In the 
angular directions, we either imagine that they are randomly distributed, in which case their 
average separation is /ib-Rb/'^b^ (the power of 1/5 is due to the flve-sphere) which is much 
larger than the local red-shifted string length along the extra dimensions, /i^m"^; or we 
imagine that they stick together and roll with a flxed angular coordinate, which is different 
from them forming a higher dimensional brane and expanding around the center. In both 
cases, we expect the leading effects of a large number of branes on density perturbations to 
be well represented by the Abelian action. As in Ref. [3S1EB], we will use this approximation 
in this paper. 

Lastly, there is a trivial slow-roll region in IR models if we tune f3 to be near (9(0.01). 
We skip this parameter space in this paper and start from /5 > (9(0.1). 

2.3 Open questions 

We list two open questions that are relevant to our parameterization: 

• Construction of potentials: Different parameter regions in the phase diagrams have 
different requirements on the inflaton mass. As mentioned, the generic magnitude of the 
inflaton mass is expected to be of order H, but the actual value can be environmental. For 
the slow-roll phases in both UV and IR models, such a mass term has to be tuned to percent 
level of the generic value; for the IR DBI model, although the magnitude of the inflaton mass 
can be of the generic order, it has to be tachyonic, or more generally the potential has to be 
repulsive for branes. For example, the conformal coupling will give a positive mass-squared 
2H^ (through the canonical inflaton dependence in the Kahler potential). It is possible that 
such a contribution gets cancelled by others from the superpotential or Kahler potential, to 
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order O.Olif^ for slow-roll models. For IR DBI models, it has to be cancelled to a negative 
value, although inflation is insensitive to the magnitude. In addition it is easy to see that, for 
IR DBI inflation to happen, the shape of the repulsive potential can be much more general 
than the quadratic form [57]. For the UV DBI model, the requirement seems to be more 
restrictive. The typical potential is expected to vary over A0 < rris ^ Mpi, so to have 
a quadratic form over Planckian size A0 needs to be functionally fine-tuned. In addition, 
using the geometric conditions (12. 7p and (12. Sp . the requirement (12. 6p implies 

m > Mpi/vOl > f m, . (2.23) 

In string constructions, the mass parameters typically arise at most of order of the string 
scale mjj, with possible suppressions from factors of the Planck mass rris/M-pi and warp 
factors. So a a ^ \/ ^ aI QsjrtA is necessary, where a a is defined before (12. 4p and is due to 
orbifolding. The construction of all these potentials is an issue under active investigation 

[SHEHEniEIlEg. 

• Background D3-charges: As we have seen, to get enough e-folds in the DBI model, 
Nb ~ 10^ is enough. But as we will see later, there are interesting parameter spaces in 
DBI inflation which require very large D3-charges for the A or B-throat to fit the magnitude 
of the density perturbations. This charge can be as large as of order 0(10^^) [33]. In IR 
models, branes are always generated in a large number after the flux-antibrane annihilation, 
and this reduces Nb to O(IO^) [35]. In a GKP-type flux compactification, such a charge 
should be cancelled by the induced negative D3-charge of the wrapped D7-branes. This 
negative charge is given by the Euler number of the corresponding fourfold in F-theory. The 
explicit examples give no more than O(IO^) [63]. So far it is not clear which of the following 
possibilities is true: in terms of the density perturbations, DBI inflation is extremely fine- 
tuned or even not viable; a modified construction, e.g. multiple-dimensional orbifolding (a 
large a^), can be engineered; a more complete understanding of the flux compactification 
can give such numbers; or subtleties are involved in the reheating. 

The approach in this paper to both issues above is phenomenological. By parameterizing 
and comparing them to experimental data, we can hopefully learn something useful about 
string theory from a bottom-up approach. 

3 IR DBI model 

In this section, we summarize the main results and predictions of the IR model, carrying 
out numerical calculations whenever is necessary. In the next section, we compare them to 
observational data, constrain microscopic parameters and make predictions. 

3.1 Attractor solutions 

In regions where various back-reactions to the warped background are negligible and the 
Hubble energy stays below the red-shifted string scale, the low-energy dynamics of the in- 



15 



flaton branes is described by the DBI-CS action 




S=Mk I d^a:,/^R- I d'xy^{^Jl+'-^g>^^d,<pd,<P-^ + V{<P)] , (3.1) 



where V{(j)) is given by (12.131) . The branes start from the tip of the throat and end at the 
UV end of the throat (pcnd = R-bV'iT'bTs = \/Xb/Rb- After that some of them quickly find 
antibranes and annihilate, diminishing the cosmological constant Vol3 

The dynamics of the inflaton can be approximately described by two attractor solutions. 
The first is that of the IR DBI inflation. This is the phase where the effect of the speed-limit 
is important. The inflaton is traveling near the warped speed of light, and the attractor 
solution is 

' ^ ■ ■ (3.2) 



t 2/3mH^ 

where t is chosen to run from —oo for convenience. Recall that for < /? ^ NB/{aBnB), H 
is approximately a constant. This phase ends around (p ~ H^/Xb {t ~ —H^^) for /5 > 1 and 
~ PHy/Xs {t ~ —f3^^H~^) for (3 <1. The inflationary e-folds as the function of can be 
estimated as 

rDBI _ H^/\^ 



^DBi ^ _ 13-^ . (3.3) 

Here we have incorporated both the case (3 < 1 and /? ^ 1 by adding the term j3~^\ this 
correction is negligible for /3 > 1, and the validity of (13.21) requires N^^^ I3~^. 

The second attractor solution describes the nonrelativistic rolling where the inflaton 
velocity stays far below the speed-limit. In this limit the equation of motion 

+ 3f/0 + d^V{<P) = (3.4) 

has the following attractor solution 

= 0oe"^^(*-*o)/3 ^ ^ -9 + V81 + 36/j _ 

Zp 

The consistency condition that the inflaton velocity is non-relativistic <C 0^/ -\/Xb requires 
that ^ aPH\/\B/3. This phase is smoothly connected to the previous DBI phase. The 
total number of inflationary e-folds provided by this period is given by 



ap 



V^B I Rb 2 



IuHRbI . (3.6) 



^ Among all the branes rolling out from the B-throat, only those which annihilate antibranes have sig- 
nificant contributions to the density perturbations. We will denote this number as ub- More generally, 
antibranes that get annihilated can reside in different A-throats. Due to different warp factors, each annihi- 
lated brane pairs can have different contributions to reheating energy. These subtleties will only affect the 
microscopic interpretation of the parameters (such as nAh\). 
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Figure 4: Attr actor solutions and numerical results. The dashed lines are the analytical attr actor 
solutions. The solid lines are numerical solutions with different initial velocities. The upper-left 
panel shows the evolution of the ratio of the inflaton- velocity r to the warped-speed-of-light h? . The 
two dashed lines are DBI and non-relativistic rolling, respectively. The upper-right panel shows 
the evolution of the Lorentz factor 7. The lower panels are the blow-ups of the upper panels. The 
parameters are /? = 2, Nb = 10^, ub = 10^, uisgs = 10~^Mpi, nAh\ = 1. In the simulation, 
branes are started at Hb = 2.9 x 10""^. 

We emphasize that this nonrelativistic rolling region is slow-roll only if /3 <^ 1, while the 
above formulae are valid even if this condition is not satisfied. For example in Fig. [3] at 
around /3 ~ 1 (m ~ y^/Mpi), after the DBI phase it takes time for branes to go through 
the lightly-shaded region till it reaches the end of the throat. This is because the branes 
are originally very close to the top of the potential. It provides a certain amount of e-folds 
typically not corresponding to the scale of the CMB. These additional non-relativistic non- 
slow-roll inflationary e-folds is also interesting to us, because it affects the relevant e-folds 
in the DBI phase, 

iVe = Nf^' + NZ , (3.7) 

and hence predictions for observations. 

In Fig. m we demonstrate numerical results and show that the two attractor solutions 
(13.21) and (13.51) give good analytical approximations for the inflaton dynamics. Any initial 
angular motions will also be damped out due to the Hubble friction, because the inflaton 
potential considered here has only radial dependence. 
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Klehanov-Strassler throat 

In this paper, we use the AdS^, geometry with a length scale R to represent the warped 
space. The details of the geometry can be different for more realistic cases. We expect our 
example to capture the main properties of the model, and to be a good approximation for a 
certain generic parameter space. Let us consider, for example, the KS throat. 



h{r) 



4^ 



12 



tot 



3^ 
27r 



In 



where we have defined a running Ri 



Rf 



with = MK,s = M{K,,, - I) for 



277r 



exp(-2/7r/35(^M) 



(3. 



(3.9) 



(3.10) 



So instead of the parameter Nb-, here we have the parameter M. The effective N^s and K^s 
are now functions of r or /. From fl3.10p we can estimate that, during IR DBI inflation, 
A/ ^ QsM. Therefore as long as 



i^eff > 9sM 



(3.11) 



we can neglect the running of A'efr. The Nb in our analyses thus represents the N^^ in a 
small region of r relevant for IR DBI inflation. The condition (13.111) is most easily satisfied 
by having a small Qs- In addition, since the WMAP window is only a few e- folds, which 
corresponds to A/ ^ ^iggM/N^^^, we only need K^^ > g^M to approximate A^^cfr as a constant 
in this window. In this case, the running of N^s only slightly affects the total DBI e-folds, 
and hence the relation in (13. 7p . Another difference between the KS throat and the geometry 
we use is that, in the latter, the UV edge of the warped space is cut off and glued to the 
bulk at i?, while in the former it is given by an independent parameter r^ax- This does not 
cause too much difference in the analyses, since the non-relativistic fast-roll inflation mostly 
happens near the top of the potential; hence A^tot^ is insensitive to the cutoff in generic cases. 



3.2 Power spectrum 

We first look at the density perturbations in the DBI phase. Its amplitude is given by the 
usual formula 



(3.12) 



where 5t is the position-dependent time delay caused by the frozen quantum fluctuations of 
the inflaton, 50 = H/2ti. In this phase we approximate the inflaton zero-mode velocity as 
the speed of light, = 0^/vXb. So we have 

(ArpBi)4 



47r2A, 



(3.13) 
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If this is responsible for ^ 23 x 10~^°, we need ~ 10^'^. Since the number of branes 
created after the flux-antibrane annihilation can be as large as 0{\/Nb), this requires Nb ^ 
10^. The formula fl3.13p can be derived rigorously using the formalism of Garriga and 
Mukhanov [6l], where we can see that the main difference from the slow- roll case is the 
development of the sound speed Cg on the world-volume of the inflaton branes. This shrinks 
the Hubble horizon by a factor of Cg. The underlying physics can be most easily understood 
in the view of an instantaneous co- moving observer with the brane [Ml [35] . For this observer 
the Hubble expansion rate is increased by the Lorentz factor 



DBI 



7 = 1/c. = 1/^1 - Ab0V0' ^ (3.14) 

due to the relativistic time dilation. In the last step of fl3.14p . we have used the IR model 
solution (13.21) and (13. 3p . This Hubble parameter leads to a horizon of size CsH~^, which lies 
orthogonal to the brane velocity and hence appears the same to the lab observer (i.e. the 
observer that does not move with the branes). 

In this model it is very important to realize the validity condition for the field theory 
analyses of density perturbations, and make estimates for the density perturbations when 
the field theory analyses break down [3ll[35l[38]- There are the following several interesting 
regions as we extend the inflaton back in time towards the IR side of the warped space. 

Firstly, open strings on the inflaton brane will be created when the Hubble energy density 
^4^4 ^Yi^Q moving observer becomes larger than the red-shifted brane tensioi]j hj^T^ = 
(^'^/ubXe- Using (13. 3p and (I3.14p . this happens at the critical e-fold 

,1/8 ^1/8 

N ^ (3 15) 

Another observation that also indicates that we cannot naively extend the field-theoretic 
results too far down the IR side of the throat, is to look at the region Ng > Nc, where the 
brane fluctuations in the transverse directions become superluminal. This is impossible. The 
reason that such a superluminal speed even occurs under the DBI action can be understood 
as follows. When we calculate the primordial fluctuations, in the first step, the source of 
such fluctuations is the uncertainty principle, which a priori does not necessarily respect 
the speed-limit if we only consider the scalar field. In the next step, the later evolution of 
such fluctuations is governed by the DBI action and always follows the causality constraint. 
The superluminal fluctuation speeds to which we just referred come from the first step, if 
we naively extend the field theoretic calculation to the regions beyond (I3.15p . 

Secondly, closed strings will be created when the Hubble energy for the lab observer 
becomes larger than the red-shifted brane tension. This happens at 

Ne ~ XT/^T ■ (3-16) 

Finally, when the closed string density created by the Hubble expansion overwhelms 
the source (fluxes or branes) for the warped geometry, the warped space gets cut off. As 



^If we replace the brane tension with the string scale m^, we have an extra factor of gl^^ in (|3.15p 
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mentioned, this back-reaction determines the maximum number of inflationary e-folds in the 
DBI phase, 

- ^T/nT ■ (3.17) 
It is important to note that the zero-mode dynamics of the inflaton are still valid as 

rDBI 

"tot 5 



long as A'^e < N^f^, since it only relies on the existence of the speed-limit and therefore on 



the condition (13.171) at which the warped space gets cut In addition, the strings and 
graviton KK-modes are only created in the tip of the throat and have energy density 0{H'^). 
It does not backreact significantly on the Hubble expansion. However, the field-theoretic 
calculation of the density perturbation is no longer valid if A'^g > Nc-, since not only scalar 
fields but also open strings will be created. 

While a rigorous treatment is currently unavailable, there are a couple of ways to esti- 
mate the density perturbations in this situation [351 [38] • (We shall make the estimates more 
quantitative in Appendix [Al) We can estimate the part of energy that goes into scalar fluc- 
tuations to be saturated when the Hubble temperature reaches the brane tension (P'^/ubXe 
at (13.151) . Further relative increase of the Hubble energy excites strings and branes. The 
stringy excitations will be diluted by the exponential spatial expansion after the Hubble en- 
ergy drops below the brane tension as branes move to the UV side, in the same way that the 
inflation dilutes relic densities. Only the scalar fluctuations are frozen and later translated 
into the position-dependent time-delay for the reheating. For the moving observer, the scalar 
field energy density is (5</))^oy7^if^ ~ (P'^/ubXe and for the lab observer 6(f) = S(f>mov/l- This 
estimate leads to 6(f) = \/XbH/ ^/nB{N^^^Y'y^ . So for A^e > A'c, we estimate 



This is also the result that we will get by looking at the transverse fluctuation speed of each 
brane. The field-theoretic analyses lead to the fluctuation speed (for the moving observer), 
^mov = 6r/6t ^ ^^/(^'yH)~^ , which is the fluctuation amplitude divided by a Hubble time. 
This velocity reaches the warped speed of light precisely around (I3.15p . Above the phase 
transition, we assume the fluctuation speed saturates the warped speed of light h"^, so Sr^ov ~ 
h?{^H)~^. The position-dependent time-delay is then 

<^iiab « ^^^^ ^ r'H-'/V^ , (3.19) 



where tq = h\ zero-mode brane speed, and the factor 1 / ^/nB is due to the reduction 

of the root-mean-square of the fluctuations by the superposition of independent branes. 
Eq. fl339|) reproduces ((SZHj). 

3.3 Regional large running of spectral index and phase transition 

From the last subsection, the spectral index is 

dhiPk 4 



^ It will be interesting to understand better how branes move through a gas of strings and graviton KK 
modes, whose effects are ignored here. 
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for N^^^ < iV„ and 



- 1 - ]^ (3-21) 

for N^^^ > N^. 

The most interesting information from (13.151) is its smallness due to the power 1/8. For 
~ 10^^, X^^^ ~ (9(100), which aheady makes interestingly smalL Considering the 
more reahstic multi-brane case ub ^ V^b leads to even smaller A^^ of order (9(10). So 
such an interesting phase transition may well have occurred within our CMB scale. Another 
interesting property is the fact that Eq. (I3.15P is independent of the inflationary energy scale 
and the local warp factor of the infiaton branes, so it will be a rather generic prediction of 
the IR DBI models0 

It is very important how sharp this transition is in terms of e-folds. In Appendix [Xj we 
give an estimate of the transition width based on the following approach. For the familiar case 
of field-theoretic density perturbations, the super-horizon perturbations can be understood 
as being generated by the random walk of the transverse brane fluctuations within a Hubble 
time before the modes exit the horizon. Such a random walk velocity is given by the Hubble 
energy, and is non-relativistic. As we have discussed, during or above the stringy phase 
transition, the main difference is that the Hubble energy is comparable to or exceeds the 
rest mass of the brane in a Hubble-size patch. As a consequence the brane fluctuation speed 
becomes relativistic. We therefore use the same physical picture underlying the familiar 
theory, but generalize it relativistically to estimate the behavior of the density perturbations 
across the phase transition. The result is given in (lA.4p . 

It is worth noting that this scenario has marked differences to several other cases com- 
monly discussed in the literature. Firstly, this model has a scale- varying running of Ug, in 
contrast to the commonly investigated empirical ansatz, where the running of Ug is assumed 
to be constant. Here the large running of Ug is only regional, principally when A^^^^^ < A'^c- 
Secondly, this scenario can generate large running, in contrast to most slow-roll scenarios. A 
standard slow-roll potential predicts very small running of Ug. A transient, large dug/dlnk 
can be caused by some "mild features" in the potential. For example, for small field slow-roll 
inflation, to generate a large transition for Ug from blue to red, the mild feature should be 
a potential shape changing from concave to convex. We study this case in Appendix [Bl 
Since the spectral index is still close to one, the slow-roll parameters for this case should 
still be at least of order 0.1. So such a case predicts unobservable non-Gaussianities. Lastly, 
non-standard choice of vacua may also cause observable running of spectral index. As we 
will discuss in Appendix [Cl such a running will be oscillatory in the WMAP window and 
phenomenologically distinguishable from the phase transition in IR DBI inflation. 



In the STA model [351 [33], when the branes move towards the IR side of the warped space, the Hubble 
energy drops linearly as (j), H k rrKjj/Mpi, the same as the warp factor. Comparing the relativistic Hubble 
energy sa jm(j>/Mpi with the red-shifted brane tension T^^'^hj^ — 0/(riA'^A)^^^, the phase transition 
happens for 7 > M-p\/ (ran'jl'^ N^^'^), in which the spectral index transitions from n,j — 1 ~ at large scales 
to — 1 ~ — 8/p at small scales. (In this footnote, we are ignoring the issue of UV embedding discussed in 
Sec.[2ll) 



21 



3.4 Large non-Gaussianity and small tensor mode 



The three-point function of the scalar perturbation for general single field inflation models, 
where the Lagrangian is an arbitrary function of X = —\g^^d^(l)dv(t) and (f), is derived in 
Ref . [65] . In the absence of any sharp features [29] , large non-Gaussianities can arise if the 
sound speed <^ 1 or another quantity A/S 1 (related to the third derivative of the 
Lagrangian with respect to X) . This non-Gaussianity is a function of three momenta, which 
are conveniently referred to as the shape of the non-Gaussianity [66] and the running of the 
non-Gaussianity [38]. The former describes its dependence on the shape of the momenta 
triangle, and the latter the overall size of the triangle. In the absence of sharp features, 
the running is relatively weak, and the shape has two categories: (1) the "local shape" in 
which the non-Gaussianity blows up in the squeezed-limit (where one of the momenta goes 
to zero) and takes a minimum value in the equilateral-limit (where all three momenta are 
equal); (2) the "equilateral shape" in which the non-Gaussianity vanishes in the squeezed- 
limit and reaches maximum value in the equilateral limit. The primordial non-Gaussianity is 
considered to be possibly observable if its estimator /tvl is > 2 or If'^i} > 10 [68l [69] . 

where the superscript "loc" refers to the local shape and "eq" the equilateral shape. 
For DBl inflation, the result becomealj 



^ -0.32c;2 . (3.22) 

For IR DBI inflation, using the relation 03.141) . we have [3B| 

^ -0.036/32(A^°^^)' . (3.23) 

The current observational bound is —256 < < 332 [67]. Comparing (12.121) and (13.231) . 
we see that the requirements of the non-Gaussianity bound on the fundamental parameters 
are quite different. Furthermore the running of non-Gaussianities for these two cases are 
opposite, as dictated by the background geometry scanned through by the rolling inflatons. 

However, we emphasize that the above results are derived in the regime where the primor- 
dial fluctuations are field-theoretic. Therefore the results can be different when the stringy 
phase transition happens. As we will see, data analysis suggests that the critical scale kc 
for that transition lies somewhere near the largest scales. For those smaller scales, it seems 
reasonable to assume that the magnitude of Eq. (I3.23P should be smoothly modified by the 
stringy corrections. So in this paper, we will use this field-theoretic approximation (I3.23P 
and the bound If^^j^l < 256. 

For DBl inflation, the scalar and tensor perturbations can be written as follows, 

where /(</>) is the background geometry and for our case / = A^/^^. So the tensor to scalar 
ratio is 



^^The papers [7D1[5S] chose an opposite sign convention of /jvl to the WMAP convention [7TJ[Tl[n7]. In 
this paper, we quote f^L in the WMAP convention. We thank Marilena LoVerde and Sarah Shandera for 
the clarification. 
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Because V is almost a constant in this model, the fact that the scalar and tensor modes have 
different horizon sizes during inflation makes no difference to (13.251) . From Sec. 13.11 we see 
that, at the scale of N^^\ (j) ^ HRly^^^/N^^\ So we get 



:iHR 



B) 



< 



riB 



ATDBi^ Ne/aB Nb/gb ' 



(3.26) 



This is of course consistent with the Lyth bound [721 HOlllS], since the r.h.s. of (I3.26P is just 
the square of A(f)/ANe in Planck units divided by (A''^^^)^, as we can see from (12. lip . (I2.17P 
and (13. 3p . To ignore the probe brane back- reactions we need 777.^ -C A^'s/a^. So 



TTS « l/(iV°^^V) < 10-^ , 

which is unobservably small. Therefore in our data analyses, we always set rxs = 0. 



(3.27) 



3.5 Constraining microscopic parameters 

In this subsection we identify the set of microscopic parameters of the model, and list self- 
consistency constraints and the observables. As discussed in Appendix [XJ we can estimate 
the power spectrum at all scales across the transition region by the following formula. 



where 



1 + 



32TT^h%T3 



ubX 



B 



327rV^ 



(3.28) 



(3.29) 



These equations relate the fundamental parameters to the observations. We will always 
choose initial position and velocity of branes so that all the observable scales are within 
the attractor solution, namely the total e- folds is larger than the minimum requirement. 
So these initial conditions will not enter the observables. The parameters Vq and (3 = 
m^/Hl determine the scale and shape of the relevant part of the potential, XbI^b char- 



acterizes the background geometry, 

0end 



and i?R tells us where to end the inflation [i.e. at 



^ ?>2tt'^ / 27 y/Xs / Rb) ■ So these four parameters determine the zero-mode evolution of 
the spacetime background and the infiaton dynamics in terms of the number of inflationary 
e-folds to the end of the inflation. In particular this determines the evolution of 7, (j)/ y/ns 
and H in (13.280 and (I3.29p . Note that we can write the factor /i^Ts in (13.290 in terms of 
and XbI^b = 2h'^Nb, h^T^ = [(j)/ ^/ubY / {\B/nB), so this is also determined. Be- 
cause the factor 0^ appears in the denominator of (I3.28p . the parameter hb affects the overall 
scale of Pfc, but does not affect the spectral index once Xb/ub is fixed. In conclusion, we have 
five parameters {Xb.ub.Rb.Vq, [3}. Using (12. 4p and (I2.18p . these parameters are equivalent 
to five even more fundamental microscopic parameters {NB,nB,msg^^^^ 
We have the following observables: 



nAh\,(3}. 
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1. The amplitude of the power spectrum ^ 23x 10~^°. Through fl3.13p and (13.181) this 
roughly determines the order of magnitude of the parameters or UbP^ depending 
on whether the pivot point iVg is smaller or larger than Nc- 

2. The scale-dependence of Pk- This determines the spectral index and its running. Since 
the spectral index of this model has a regional large running, the data will constrain 
at which scale (kc) such a running happens and which DBI e-fold (Nc) it corresponds 
to. These are then transferred into some delicate relations with the microscopic pa- 
rameters, e.g. (1A.7P and (13. 7p . 

3. Non-Gaussianity constraint 7 < 28, will mainly constrain [3 (with some weak depen- 
dence on iV°^^). 

4. We have the following several consistency relations. First, a scale k is related to the 
corresponding by 

iV, = 65 - In + In , (3.30) 

0.002 MpC^l Treheat 

where the reheating is assumed to be efficient!^ so that Trchcat = Vq^'^, and is the 
sound speed when the mode k = 0.002 Mpc~^ crossed the sound horizon. 

Second, according to the multi-throat brane inflationary scenario, the maximum num- 
ber of the infiaton branes is bounded by the flux number M. Since Nb = QbKM and 
we want to keep the minimum warp factor (I2.22p small, we require 



riB < VNB/{aB9s) • (3.31) 



Third, the geometric constraints (I2.7P and (12.80 give an upper bound on nisg. 



-1/4 



^<2='V'M^^. (3,32) 

where the approximate numerical factors come from the toroidal compactification, 
Mpj = j^^^g~'^m^gV^, which may change for more realistic setups. 

Fourth, the warp factor < 1 and ua = ub, (this is not a coincidence, see footnote 

CD, so 

nAh\ < ub ■ (3.33) 



Lastly, the inflationary scale and the string scale are both bounded below by TeV, 

UAh^Ts > TeV^ , (3.34) 
rrisgj^^^ > TeV . (3.35) 

^^For single-throat reheating, the brane-antibrane pairs immediately (in terms of the Hubble time) anni- 
hilate and decay into relativistic particles and start the usual radiation-domination epoch. For tunneling 
reheating such as double-throat reheating, (|3.30p may receive some small modifications due to a long inter- 
mediate matter-domination epoch |51j . 
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(These two bounds turn out to be very weak. Much stronger ones will arise from the 
data analyses.) 

Note that, there are two other independent parameters Qs and as that only appear 
in the bounds (13.311) and (I3.32p . For simplicity we do not promote them into free 
parameters in data analyses. We set Qs = Q.l and = 1 in these bounds. Reducing Qs 
and/or increasing as may loosen these bounds and allow some microscopic parameters 
to take wider ranges. But this should not change the model predictions qualitatively. 

4 Markov Chain Monte Carlo data analysis 
4.1 Methodology 

In order to obtain multi-dimensional parameter constraints from cosmological data, a Markov 
Chain Monte Carlo (MCMC) approach [73l[71l[ISl[l6l[771[78l[79]is employed to sample the 
likelihood surface efficiently. The MCMC is used to simulate observations from the posterior 
distribution P{a\x), for a set of parameters {a} given an event {x} (which, for us, is the 
total set of observational data), using Bayes' Theorem 



where P{x\a) is the likelihood of the event x given the model parameters a, and -P(a) is 
the prior probability distribution of obtaining a model parameter value a. The MCMC 
generates random draws from the posterior distribution that are a "fair" sample of the 
likelihood surface, and from this sample, we can estimate all the quantities of interest about 
the posterior distribution (mean, variance, confidence levels). 

In most cosmological analyses, fiat priors, P(a)=constant, are assumed on a set of em- 
pirical parameters such as the spectral index and its running, n^, dus/dlnk, and the nor- 
malization As, of the primordial scalar power spectrum, or its logarithm InAg (for example 
[27| [80l [H [m [25]). It is by no means true, however, that such constant priors should natu- 
rally arise in a fundamental theory. The effect of priors on constraints on slow-roll infiation 
was recently discussed in here we discuss their role in IR DBI infiationary scenarios. 

Unlike parameters used in an empirical ansatz, the relationships between the fundamen- 
tal microscopic parameters and observables are highly nonlinear and far from transparent. 
This can make it problematic for the MCMC to efficiently explore the likelihood surface, 
potentially leading to the presence of non-Gaussian posterior distributions: for example, 
multiple, disconnected maxima in the likelihood surface, or long, curved degeneracy direc- 
tions. In these proposal distribution for the microscopic parameters that samples the 
posterior distribution efficiently can be very difficult to obtain. In such situations, instead 
of directly adopting the microscopic parameters as the parameters sampled by the MCMC, 
we find it useful to reparameterize variables according to the properties of the models. The 
specific details will of course be model-dependent, but there are certain general strategies 
that one can follow, which we will summarize now. 




(4.36) 
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• Although the full relationships between the observables and microscopic parameters 
{a} are usually complicated, and in realistic cases often have to be computed numeri- 
cally, an isolated analytical expression for the observationally accessible window (scales 
10-^ Mpc"^ <k<l Mpc-^) can be much easier to obtain and be expressed in terms 
of an equal or smaller number of effective parameters {6}. 

• Run a trial MCMC with the effective parameters {9} with constant priors, in order 
to ensure that these parameters have a relatively simple likelihood surface. This will 
generally be the case if the {6} are chosen to such that the observables of the model 
vary roughly linearly with the effective parameters. 

• The effective parameters {6} can often provide the necessary physical intuition to find a 
reparameterization of the original microphysical parameters ai{ai), which have simple 
(e.g. linear) relationships to the {6}, thus have simple enough relationships to the 
observables such that the likelihood surface can be effectively explored by standard, 
robust MCMC techniques. Ideally, the reparameterization Q!j(aj) should be a bijective 
function in the observable region of interest. Because the trial MCMC helps ensure the 
simplicity of the likelihood surface in the space of {6}, and the new parameterization 
ensures that the {a} essentially travel along the directions of {6}, the likelihood surface 
in space of {a} will also be plausibly simple. When running the full MCMC in the 
{a} space, any analytical approximations used in the trial MCMC to compute the 
observables can be dropped, and the observables calculated numerically, in order to 
prevent modeling uncertainties coming from such approximations significantly affecting 
the final constraints. 

• After obtaining the likelihood surface of the new parameters {a}, transform the like- 
lihood surface of the {a} to the space of the original parameters {a}; the MCMC can 
also be re-weighted to impose any desired priors on the {a}. It must be noted at this 
stage that the theory does not predict the prior distribution of the {a} and therefore 
any prior adopted on this parameter set can potentially be highly informative. If the 
data impose a tight-constraint on a given parameter {i.e. the likelihood is significantly 
peaked within the prior) such that the posterior distribution is not very sensitive to 
simple forms of adopted prior (such as constant, logarithmic etc), we will not concern 
ourselves overly with this point. If a given "constraint" is coming primarily from the 
prior, we will point it out. 

• An alternative approach, which is to use complicated sampling techniques to explore 
the complex likelihood surface of the original microphysical parameters {a}, can often 
be more time-consuming as it has to be tuned for each particular problem. 

Now we will apply this procedure to the IR DBI model. 
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4.2 MCMC using microscopic parameters of the IR DBI model 

First, as detailed in Appendix \^ we notice that the primordial power spectrum predicted 
by the IR DBI model can be approximated by the following analytical form: 



R = H'Sf '''''' 



1 - 



N, 



16 



(4.37) 



This is parameterized by three effective parameters: Nc, In/c^; and nsP'^, where kc is the 
critical scale near which the stringy phase transition happens, 

iVf = ln(A;e/A;) + A^, . (4.38) 

After verifying that these three parameters appear to have a simple likelihood surface in 
a trial MCMC analysis, we relate the five microscopic parameters to these three parameters 
through approximate analytical expressions. 

The relation to A^ is simple and given by (lA.TD . 

= Ven^'^ ■ (4.39) 

The kc is defined as the value of k at A'^^^ = A'^. Using the relations (13. 7p and (13.301) . 

we get 



0.002 Mpc'V VCsTreheat 

Using the approximation (13. 6p . expressing Hq, Rb, Vq in terms of A^^, n^/i^ and gs/^t using 
Eq. ([MD and flZTSj) . we obtain 

k \ 13 a^^'^ 3 

^ ^ 65 - A'c + In \ In \ m imb 

0.002 Mpc^V V67r3/4c, a(3 ^t:^!^ Aaf3 



ln( -] ^ 65- Nc + ln^^ + — In — + — ^lnA^£ 

VQit'^/^Cs a(3 y6vr5/4 4a/5 

l,_i_),„„,^._(l,_i_),„^, (4.41) 



s 



where 



a(3 = (-9 + v^Sl + 36/3)/2 . (4.42) 



Here, is the sound speed when the mode k = 0.002 Mpc^^ crosses the sound horizon; 
we can also approximately express it in terms of the five microscopic parameters. But the 
detailed expression will complicate the relation. For our purpose, since Cg varies slowly from 
0.01 to 0.1, treating it as a constant should not cause problems for the reparameterization. 

These expressions suggest that the following set may prove to be a successful reparame- 
terization of the microphysical parameters that can be effectively explored by MCMC: 
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as = InnsP^ . (4.43) 

The relation between these new parameters, {a}, and the effective parameters, {6} = 
{Nc,lnkc,nB(3'^}, is very clear. Two of them are identical, and the rest, a2, a^, and a^, 
all have approximately linear relationships to In kc through (14.411) . 

We adopt the reparametrized microscopic parameters {a} and the standard set of cos- 
mological parameters {ub = Vt^h?, uom = ^mh^-, QaiT} as the model parameter set sampled 
by the MCMC. Here, 9a is the angular size of the acoustic horizon and functions as a proxy 
for the Hubble constant Hq = lOO/i km/s/Mpc or Qm, and r is the optical depth to reion- 
ization. The universe is assumed to be spatially flat. Constant priors are assumed over the 
previously specified parameter set {a}, subject to the microphysical cuts described below. 

For each set of {a} sampled by the MCMC, the relations (14.431) are numerically inverted 
to obtain the set of microscopic parameters {Nb, ub, nAh\, Qs ^^'^rrig, (3}. This inversion is 
bijective in the parameter ranges of interest. The microscopic parameters are then fed into 
a numerical code which is described in detail in Appendix [Dl After checking the input pa- 
rameters for a set of microphysical conditions which enforce model-building self-consistency, 
as described in the Appendix [Dl the code computes the primordial power spectrum of the 
curvature perturbation for the model specified by the input parameters. Input parameters 
which fail to satisfy the microphysical cuts are rejected through being assigned zero likelihood 
in the MCMC. The primordial power spectrum from this code is fed to the Boltzmann code 
CAMB [76], without significantly increasing the computational time, in order to calculate 
the cosmological observables. 

We use a modified version of the CosmoMC code [77] to determine constraints placed on 
this parameter space by the WMAP three- year cosmic microwave background data [1], [83l 
[m [85] and the SDSS Luminous Red Galaxy (LRG) galaxy power spectrum data [1]. We 
marginalize analytically over the linear bias factor b and the non-linearity parameter Qni of 
the SDSS LRG data as is done normally in the CosmoMC code. A properly derived and 
implemented MCMC draws from the joint posterior density defined in (I4.36P once it has 
converged to the stationary distribution. We use eight Markov chains and a conservative 
Gelman-Rubin convergence criterion |86] on the eigenvalues of the parameter covariance 
matrix to determine when the chains have converged to the stationary distribution. Then 
we re-weight the MCMC to switch to constant priors on the microscopic parameters {a} = 

{logio ^B, logio Ub, (3, logio ^AhA, logio ^s/gy^}- 

Following this process, we would like to apply the observational non-Gaussianity bound 
—256 < < 332 (95% CL) [HZ], as this should have a significant effect on restricting the 
allowed parameter range for (3 and other parameters which are correlated with it. This is 
because 7 is roughly proportional to /?, and f^j^ = —0.327^. However, two approximations 
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Model 


Best fit -21nA^ax (WMAP+SDSS LRG) 


ACDM 


5374.04 


IR DBI 


5373.11 



Table 2: The best fit chi square, defined as = —2 In/^max (where £max is the maximum likelihood 
with respect to the WMAP 3 year data and the SDSS LRG galaxy power spectrum data) for the 
standard ACDM model and the IR DBI scenario analyzed in this work. A Gaussian prior on 
has been applied based on the WMAP 3 year constraint on this parameter from Ref. [67] . For the 
ACDM model, the /atx, constraint has been applied assuming = 0. The IR DBI model gives a 
slightly better (lower) foi' this dataset than the ACDM model. The primordial power spectrum 
is described by five microphysical parameters in the former, and two empirical parameters (an 
amplitude and a power law index) in the latter. When we consider that the IR DBI observables 
are described phenomenologically by the three effective parameters Nc, InA;^, and n^/?^, to which 
the microphysical parameters are related, we can see that the IR DBI model has roughly one extra 
degree of freedom over the ACDM model, which one expects to give a Ax^ ~ 1 improvement in 
the fit. Since this is in fact what we see, there is no indication of a preference in the data for the 
IR DBI model. 

enter in applying this constraint. First, the observational constraint was obtained using an 
estimator that does not encode the specific scale-dependence of the IR DBI model, and it also 
does not restrict < 0. Second, Ref. [67] only gives a 95% confidence level of the result, 
and hence a full Bayesian posterior is not available for this parameter. In order to make use 
of this constraint despite these limitations, firstly we assume that the constraint of Ref. [67] 
is the effective constraint at k = 0.02 Mpc~^, which is approximately the best constrained 
scale with the current data compilation [871 EH]- Secondly, we choose a Gaussian prior 
on f'^j^{0.02 Mpc~^) which has the 95% CL range found by Ref. [67], since the maximally 
uninformative prior in the case that only a single confidence range is available has a Gaussian 
form [SHI [90]. We apply this prior to the chains, verifying that the convergence criteria still 
remain satisfied. Finally, we obtain parameter constraints on the microphysical parameters, 
cosmological parameters, cosmological observables, and derived model parameters, which we 
present in Tables |2H3] and Figs. [514TQ1 
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Parameter 


Marginalized Constraint 


Maximum Likelihood 




n 091 4c;+*'-''*'071+0. 00138 


0.02162 




U.iU (U_o.0044-0.0082 


0.1058 


T 


U.UOJ_o.030-0.061 


0.094 


Ho 


'71 9+1.8—3.9 
' J^-^-1. 9-3.7 


72.1 




A «/| +0.30+0.45 
^•"^-0.32-0.70 


4.93 


174 

^ogioirris/gs' /Mpi] 


71 +1.04+1.43 
^- ' J--l.07-2.89 


-5.91 




9 1 1 +0.63+1.63 
^•-'--'--0.60-0.85 


1.77 


loglo[iVB] 


Q /IO+0.39+0.93 
^•^°-0.39-0.70 


9.15 




1 1 +2.64+3.34 
J- •^-'--2.82-8.02 


0.585 


logio[K'/VMpi] 


nc;+i-25+i.83 

^•-"^-1.34-3.09 


-6.36 




qcr 7+6.8+11.7 
' -7.3-12.6 


34.1 


logio [kc/Mpc] 


/I 1 e;+0.81+1.21 
^•-■-^-0.81-1.82 


-3.86 






20.5 


N^^^ (10-VMpc) 


00 /1+5. 6+9.1 
'-"-'•^-6.1-10.7 


37.6 




> -2.36 (95% CL) 


-1.09 


log,,[ihAmJgyy/il6n^)/M^,] 


1 7 9+2.4+3.5 

' -^-2.6-6.0 


-16.2 


n, (0.02/Mpc) 


n Q/1Q+0. 016+0.032 
U.a^O_Q 016-0.031 


0.946 


dus/dlnk (0.02/Mpc) 


n noi +0-008+0.011 

U.UZl_o.009-0.025 


-0.021 


7 (0.02/Mpc) 


1 Q Q+3.6+9.3 
-'--'•-'-3.4-5.1 


16.8 


(0.02/Mpc) 


1 01 +44+61 
-L -J -"--45-141 


-91 



Table 3: Constraints on the IR DBI model from the WMAP and SDSS LRG data-sets (mean, 
upper and lower 68% and 95% CL, marginalizing over all other parameters), and the maximum 
likelihood values of the parameters found in the MCMC. A Gaussian prior on has been applied 
based on the WMAP 3 year constraint on this parameter from Ref. |67j. 
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Figure 5: Solid lines show the marginalized 2D-joint 68% and 95% probability contours (off- 
diagonal panels) and ID marginalized probability distribution (diagonal panels) for the micro- 
physical IR DBI parameters. The color coding in the off-diagonal panels shows the marginalized 
probability density in these 2D parameter spaces, ranging from red for the highest density to blue 
for the lowest. 
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15 ^2,0 25 30^^, 40 50 -4 -2, 

Nf (10-^/Mpc) 10910 [(V>b)^''] 

Figure 6: Marginalized posterior probability distribution functions obtained from the MCMC 
analysis for observables and derived quantities of interest. The functions are normalized such that 
the area under the curve is one. 
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Figure 7: Examples of 2D contours. Solid lines show the marginalized 2D-joint 68% and 95% 
probability contours for observables and derived quantities of interest. The color coding shows the 
marginalized probability density in these 2D parameter spaces, ranging from red for the highest 
density to blue for the lowest. 
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WMAP 
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Figure 8: Reconstructed 68% (dark) and 95% (light) CL constraints on the primordial scalar 
power spectrum for the IR DBI model. The range of scales spanned by WMAP and SDSS LRG 
data (which were used in the fit) and the smaller scales covered by Lyman-a data (which were not) 
are shown for reference. For comparison, the dashed lines show the corresponding 68% and 95% 
constraints for (Upper) single-field slow-roll inflation, taken from Ref. [87] fitted to the WMAP 
and SDSS main galaxy sample data [3] , and for (Lower) the empirical power law ansatz where the 
primordial power spectrum is described by its amplitude at a pivot scale, the spectral index n^, 
and its running dus/dlnk, fitted to WMAP, SDSS LRG and Supernova Legacy Survey [6] data. 
See text for discussion. 
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Figure 9: Reconstructed 68% (dark) and 95% (light) CL constraints on the scale-dependence 
of the running of the spectral index, dng/dlnk, showing a mild indication for a "running of the 
running". See text for discussion. 




Figure 10: Reconstructed 68% (dark) and 95% (light) CL constraints and mean (dashed) for the 
non-linearity parameter for the IR DBI model, as a function of scale. Note that DBI inflation 
predicts < 0, and the absolute value is plotted for convenience. The range of scales spanned 
by WMAP data are shown, as well as the WMAP 95% lower limit > —256 (solid line), from 
an analysis which treated f'^^ as scale- independent [67|. See text for discussion. 
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5 Discussion and Conclusions 

We conclude by highlighting the main results and discussing their physical implications. 
The quoted ranges are at the 95% confidence level, and we have combined constraints from 
both the power spectrum and non-Gaussianity. The detailed 68% and 95% CL marginalized 
constraints and the maximum likelihood values are listed in Table [3l 

5.1 Microscopic parameters 

• Shape of the inflaton brane moduli potential: 1.3 < (3 < 3.7. The lower bound is 
due to constraints from the power spectrum, while the upper bound is due to the 
non-Gaussianity constraint. It is encouraging that, while IR DBI infiation can happen 
for a range of (3 that varies over nearly 10 orders of magnitude, 0.1 < /5 < 10^ (see 
Eq. (12.211) ). comparison with data picks out a very small range around 0{1) which is 
generically expected theoretically. This makes an explicit construction of such poten- 
tials a more interesting question. 

• Fundamental string scale: —9.6 < \ogiQ{ms /Mpi)/^y^ < -5.3. The upper bound on 
the string scale is due to the large charge, and hence length scale, of the B-throat 
required to fit the amplitude of the density perturbations. The lower bound is due 
to the fact that a smaller string scale tends to increase the total number of e-folds 
of non-relativistic fast-roll infiation, and make the running of the spectral index too 
large (Fig. [T]). The model prefers an intermediate fundamental string scale, 10^ GeV < 
i^s/gy^ < 10^^ GeV, and therefore an intermediate large volume compactification, 
8.9 X 10'' < V^^^Mpi < 4.8 X 10^^, where V is the compactification volume. 

• B-throat charge: 8.8 < log^Q Nb < 10.4; Number of inflaton branes: 3.9 < logi^UB < 
5.1. In terms of the GKP-type warped compactification, this implies fiux numbers K ~ 
M ~ \/Nb ~ (9(10^). Explicit construction remains an open question as discussed 
in Sec. 12.31 In the multi-throat brane infiation scenario, infiaton branes are generated 
from fiux-antibrane annihilation. The number of branes generated in this process is 
roughly determined by the fiux number M. Indeed, a small number of infiaton branes 
is ruled out by the data. 

• A-throat minimum warp factor: —2.4 < logj^Q /i^ ^ 0. This is from combining the 
constraint on and n^/i^, = {uAh^/nBY^^. A smaller Ha leads to larger A^^^^ 
and larger running of the spectral index (Fig. [7]). So the A-throat tends to be short. 
This makes tunneling reheating possible, where many interesting phenomena can occur, 
such as an intermediate matter-dominated epoch. 

5.2 Secondary derived parameters 

• Inflationary phases. In this model, not all e-folds comes from IR DBI infiation. The 
last 13 < N^^^f < 24 e-folds come from non-relativistic fast-rolling infiation, which is 
possible because infiatons are close to the top of the potential. 
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• The stringy phase transition. The Hubble-expansion induced stringy phase transition 
happens at the largest scales in the sky, —6.0 < log^Q kc/Mpc < —2.9. However its 
impact on density perturbations extends over to shorter scales, such as generating a 
transient large running of the spectral index. 

• Inflation scale: —10.0 < logi^V^^^ /Mp\ < —5.1. This gives a very small tensor to 
scalar ratio rrs < 10"^^. 

1/2 

• Cosmic string tension: —23 < log^gGyUD + logiofi's < —14. Here the cosmic strings 
refer to the D-strings left over from the brane-antibrane annihilation in the A-throat, 
whose tension is G/id = {mshA/ gy'^MpiY /{IGir'^gy^). There is an unconstrained free- 

1/2 

dom coming from the additive factor \ogiQ gs , but it is not expected to give any 
significant contributions. The F-string tension differs by a factor of gs, fir = Osf^D- 

5.3 Observational predictions 

• Large, but regional, running of spectral index: —0.046 < dns/d\nk{k = 0.02/Mpc) < 
—0.010. A reconstructed full-scale power spectrum and the running of the spectral 
index are shown in Fig. [8] & [91 

This prediction is stringy in nature. A better understanding of the theoretical details 
and better measurements of both the power spectrum and non-Gaussianities on the 
relevant scales may reveal finer structures. In future experiments, Planck is expected 
to achieve a{dns/d\nk) = 0.005 [8]. 

• Large non-Gaussianities: —272 < f^^ik = 0.02/Mpc) < —70. A reconstructed full- 
scale prediction is in Fig. [101 which shows the running of the non-Gaussianities. 

This prediction is strictly speaking field-theoretic, but with strong string theory moti- 
vations, such as warped compactification and the DBI brane action. This field theoretic 
regime is k > k^, the theoretical analysis for non-Gaussianities at < is currently 
unavailable and remains an interesting open question. In future experiments, on CMB 
scales, Planck can achieve cr(/^^) = 67 [HH [68]; on large scale structure scales, some 
high-z galaxy surveys can reach similar or better precision [69] . 

As seen from these results, constraints from cosmological data, and even relatively loose 
constraints such as the non-Gaussianity constraint, are already putting strong restrictions on 
models which aim to provide self-consistent microphysical descriptions of the early universe. 
With the bounty of precision cosmological data expected in the future, the hope of probing 
not just field-theoretic, but string-theoretic early universe physics burns brightly. 
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A Estimate the effect of phase transition on spectral index 

In this appendix we estimate the transition behavior of the spectral index between two 
asymptotic values described in Ref. [35], [38] and Sec. 13. 2[ Consider a simple model where the 
density perturbations are caused by the scalar field fluctuations, which are the super-horizon 
ripples on branes in transverse directions. These ripples are generated during a Hubble time 
while they are still sub-horizon and then frozen. The amplitude of the ripples are given 
by the fluctuation speed of a Hubble-sized patch on the brane. This speed is determined 
by the energy pumped into the branes by the Hubble expansion. This model simplifies 
the underlying physics by focusing on only the overall fluctuation speed of a Hubble-sized 
patch while ignoring the detailed world-volume theory such as effects from specific stringy 
excitations. 

According to the special relativity, an object with rest mass mo and energy E = mo + AE 
has velocity 



where h^T^ is the red-shifted brane tension. The Hubble energy is 'jH/2'it, half of which 
goes to the kinetic energy of the transverse oscillation of the brane AE = '-/H/ An, while the 
other half goes to the tension of oscillations in terms of spatial derivatives. We have restored 
the factor of 2-7r in the Hubble length and energy in order to quantitatively match the known 
results in the low energy limit. The local speed of light is c = h"^. The position-dependent 
time delay is 




(A.l) 



For the on-brane observer, a Hubble-sized patch has rest mass 




(A.2) 



6t 



v{-iH/2n)-^/-i 



(A.3) 
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where the numerator is the fluctuation amphtude within a Hubble time viewed from the lab 
observer (hence an extra factor of I/7 due to Lorentz contraction), and the denominator is 
the overall brane velocity which is approximately the local speed of light r ^ c. Here we also 
consider the case oi ub multiple branes where the superposition of independent fluctuations 



reduces the time-delay by a factor of 1/ 
spectrum 



Ub- Using these estimates we obtain the power 



(A.4) 
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where 



(A.5) 



This formula recovers the usual field theory result in the limit of non-relativistic fluctuation 
speed. This includes non-relativistic- (slow or fast)-roll inflation, and DBI inflation below the 
phase transition. This formula also gives an estimate on the effect of the Hubble-expansion 
induced stringy phase transition. The estimate is expected to provide the envelope behavior 
beyond the transition since it ignores detailed features such as specific resonant production 
of various stringy states. 

It will be useful to extract the DBI inflation region in flA.4p and flA.Sp . and parametrize 
it in the following way0 



3247r2 



N, 



16 



(iV8 + iVpBl8)2 



Nc is deflned as 



X 



1/8 

B 



(A.6) 



(A.7) 



where we have used the relation (12.31) . Taking the limits A'^g <^ Nc and ^ N^., we recover 
fl3.13p and f l3.18p respectively. The spectral index is 



x^ + 3x-2 



N- 



DBI 



This formula interpolates between two asymptotic values 4/A''e and —A/N^. If we deflne the 
width of the transition region as the e-fold difference AN^ between rig — 1 = 2/N^^^ and 
_2/jYDBi^ then we have 

ANe ^ 0.2^", , (A.9) 

which can be quite large (for example, six if Nc = 30). But the running of Ug is still 
observably large in the transition region (for example, dus/dlnk ^ —0.02 in the range of 
f lCTl for Nc = 30). 



More precisely, because of the sound horizon is time-dependent, we should replace N^^^ in (|A.6P with 



N^^^ — \ii{cs{k)H{k) ^)/{csH ^), where the variables with a hat are evaluated when the reference mode k 
(e.g. k = 0.002/ Mpc as in (|3.30p ) crosses the sound horizon. Because the relevant scales for (|A.6p span only 
a few e-folds, the change of the sound horizon CsH~^ is very small and we neglect such corrections. 
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B Running spectral index from slow-roll potential with mild fea- 
tures 



Usual slow-roll inflation gives negligible running of the spectral index, dns/d\\i k = O{10/N^), 
because large running of ris tends to end the inflation too quickly. For a comparison 
with data, see Ref. [92]. In this appendix, we study the possibility of having measurable 
\dns/dlnk\ > 0.01 by adding some mild features to the slow-roll potential, and how we can 
phenomenologically distinguish it from the IR DBI inflation model. 

We consider a potential of a small field inflation and add some ripples on it, 

V = Vo-a0- 6sin(0/0o) • (B.l) 

The inflaton starts for example at 0j = 0. At = (pond, one imagines that the inflationary 
energy Vq gets annihilated as in brane-anti-brane inflation. This is just an example of many 
possibilities, which we use to illustrate the main properties. As we will see, to generate 
a large running of Ug from blue to red, the shape of the slow-roll potential changes from 
convex to concave. The oscillatory ripples help to sustain inflation, and at the same time 
generate large dus/dlnk periodically. In fact, for our purpose, it is not necessary to make 
the mild feature periodic, for example the 3rd term in fIB.ll) can be regional as long as it 
falls into the WMAP range. Nonetheless being periodic might be more naturally realized in 
model-building. 

We want the inflaton to continuously roll down, so we need V < 0, i.e. , 

- a + 6/00 < . (B.2) 
We require the average slow-roll parameter 

- t i? « ^ (B.3) 

to have enough inflationary e-folds. To have the effect of one ripple span several e-folds, we 
need 0o/0 = ^H~^, where ^ is of order one. can be estimated using the attractor behavior 
3H(f) + V' = and taking the average value of V . We want the other slow-roll parameter rj 
to vary between order ±0.1 to generate observable dus/dlnk, 

where C is 0(0.1). 

Therefore for our purpose, we can choose the parameters in (IB.ip in the following way. 
We require 



so that e <^ 1 



«^ (B.6) 
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Figure 11: Turning on the periodic mild features ( does not significantly affect the total inflationary 
e-folds. In this figure, ^ = 4, a/Vo = 10~^, starts from and ends at 0.006Mpi. 



1.3| I 0.04 
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Figure 12: The power spectrum and spectral index in the WMAP range for a slow-roll potential 
with mild features. The parameters u/Vq = W~^, = 4 and C = 0.15 are chosen so that these 
observables look close to what we obtained in the IR DBI model. 
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with ^ of (9(1) so that the effect of one period of the mild feature spans a reasonable amount 
of e-folds; 



b=ec 



a 




(B.7) 



with ( of 0(0.1) so that the running of Ug is observably large. In order for these conditions 
to be consistent with (lB.2p . we need 



For example, we can choose ^ = 4, ^ = 0.15. 

In Fig. [11] and [12] we demonstrate numerically: 

• Introducing such mild features does not significantly affect, and in fact can slightly 
increase, the total number of inflationary e-folds. 

• The power spectrum within the WMAP scales looks like what we obtained in IR DBI 
model; notably, the spectral index runs from blue to red with a large and negative 
dus/dlnk. 

To experimentally distinguish this case from the IR DBI model, we estimate the non- 
Gaussianity. In slow-roll inflation, the 3-point function of the gauge-invariant scalar pertur- 
bation receives contributions from the following sources [701 [93] . In the cubic action there 
are terms proportional to e^, and edrj/dt. In terms of order-of- magnitude, these terms 
contribute C?(e), C(e^), 0(A?7), respectively, to the non-Gaussianity estimator f^^. The A?7 
is the maximum change of rj caused by the features, since the 3-point function involves an 
integration over time. In addition, there is a field redefinition term that contributes 0(?7cnd) 
to /atl, where ?7end is the frozen value of rj after the horizon crossing. 

In case of slow-roll inflation with smooth potential, the leading terms of f^L are 0{e) and 
0{ri) [To]. In case of sharp features, (9(A?7) term dominates [29]. In the case of interest here 
with periodic mild features, both 0{Ari) and the boundary term 0(?7end) become important. 
As we saw, to generate large but reasonable running of the spectral index, we require rj to 
vary between ±(9(0.1). So we expect such features to be associated with non-Gaussianities 
fNL = (9(0.1). This is clearly observationally distinguishable from the IR DBI inflation case. 

C Running spectral index from non-Bunch-Davies vacuum 

In this Appendix, we study how the running spectral index arising from the non-Bunch- 
Davies vacuum case can be phenomenologically distinguished from the running spectral 
index in IR DBI model. 

In the field theory of density perturbations, the Bunch-Davies (BD) vacuum is the leading 
behavior of the fluctuations when they are well within the horizon. However corrections to 
such a vacuum can have observational effects and may provide information on new physics 
|94j . For reviews and references see [H51 [HSj- Denoting the scale of the new physics as 
M ^ H, this correction typically arises at the order H/M, for example, if we choose the 



(B.8) 
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adiabatic vacuum at the scale M [HZ]. This is also called the trans-Planckian effect if M is 
regarded as M-p\. In this appendix, we will treat M to be much more general. 

The main difference between the case of non-BD vacuum and the Hubble-expansion 
induced stringy phase transition in DBI inflation is that, in the former, there is a large 
region between the new physics scale and the Hubble horizon where the conventional field 
theoretic analyses still applies; while in the latter, the Hubble horizon is comparable to or 
smaller than the new stringy length scale in the infiaton sector. 

The effect of non-Bunch-Davies vacuum in slow-roll inflation typically results in an os- 
cillatory modulation on the usual power spectrum [981 EH]- So it may also introduce an 
observable running spectral index. The potential for observing these features are discussed 
in Ref. |1001I101] . To study both slow-roll and DBI inflation, here we generalize the analyses 
of Ref. [971 1102] to the case with arbitrary sound speed. 

The quadratic action for the gauge-invariant scalar perturbation C, in general single field 
inflation is 



did X 



(C.l) 



where a is the scale factor, Cg is the sound speed, and e = —H/H"^ is one of the slow variation 
parameters. Using the variable = -s^Ck (-2 = o.\/2e/cs) and its conjugate momentum 



TTk = t'k - — 



the Hamiltonian in Fourier space is 



{2n) 



d^k- 



TTkVTk + k C^VkV^ + — (vTkt'k + VTi^fkJ 



(C.2) 



(C.3) 



where the prime denotes the derivative with respective to the conformal time r. 

We can quantize f k and TTk in terms of the creation and annihilation operators which are 
either time-dependent 

1 



fk 



TTk 



aklrl 



—'I 



kCa 



ak(r) 



alkl^. 



(C.4) 



or time-independent 



fk 

TTk 



/fc(r)ak(ro) + A!(r)at_k(ro) 



—I 



gk{r)ak{To) - gl{T)ali^{ro) 



(C.5) 



where //c(t) is the solution of the equation of motion v'^ + c^k'^V]^ — {z" / z)v\i = 0, 



Mr) 



9k{r) 



\/2csk 



kc 



kCgT 



^J2csk 



kc.T 



■ikCsT 



kCa 



ikcsT 



(C.6) 
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Equations (lC.4p and (IC.Sp are related by the Bogolubov transformation 

ak(-r) = 6(r)ak(ro) + Cfc(i")aU(7"o) , 



a 



= e^(r)aU(^o) + C(r)ak(ro) , (C.7) 



where 



The following relation should be satisfied to preserve the commutation relation for ( 1C.7|) . 

\^k\' - \Ck\' = I- (C.9) 

Therefore we have fkgl + f^gk = I and |C+p - |C_p = 1. 
An adiabatic vacuum |0, tq) can be chosen as 

ak(ro)|0,ro) = , i.e. ^(^0) = . (C.IO) 

The Bunch-Davies vacuum corresponds to sending tq to — 00. More generally one can choose 
a finite tq ~ —1/aoHo for the mode k, when this mode crosses the scale of the new physics 
M ^ H/cg. Hence the relation between the power spectrum in the non-BD and BD vacuum 
is 

pnonBD _ ^ C_\^P^^ , (C.U) 

where 

, ^ ^ ,., , aH\ ( 2kcs 

C+ + C\' ^ 1 - — sin 



kcs / V (iH 
1 - -— sin — , (C.12) 



where the new physics scale M = {k/a)r. The extra contribution of the non-BD vacuum to 
the spectral index is 

An, ^ -2(e + s + /i)cos^^ , (C.13) 

H 

where fi = M/{HM), s = Cs/{Hcs). In the following discussion, we will concentrate on the 
amplitude and oscillation frequency of these features on the power spectrum and spectral 
index. We want to compare them to those in IR DBI model, where the spectral index Ug — l 
changes between ±(9(0.1) within (9(10) e-folds without oscillations. For this purpose it is 
useful to note that the change of the arguments in the trigonometric functions in fIC.llI) and 
flC.lSP as a function of k can be written as 

44 





11 10 



Figure 13: Illustration of the effect of the non-BD vacuum and the stringy phase transition in IR 
DBI model, combining ()A.6P and (jC.lip . We use Nc ~ 40, which is also the scale where M ~ H/cs- 
The ripples at small scales are due to the non-BD effects and the suppression in large scales is due 
to the phase transition. The connection between them is smoothed out by hand due to a lack of 
more detailed understanding. The right panel is a blow-up of the left. 



We first look at slow-roll inflation, where s = fi = and e < 0.01. So the variation of 
Us, typically smaller than (9(0.01), is much smaller than that caused by the phase transition 
in IR DBI case, although the oscillatory frequency of the Ug is adjustable depending on the 
values of {2M/H)t-o and e. For Us to have larger variations, one needs the special case of 
e ^ 0.05; at the same time, for the running to span (9(10) e-folds without oscillation, from 
f lOTil) . we see that McjH = 0{7r). This barely satisfies McjH > 1. We conclude that 
in slow-roll inflation the effect of the non-DB vacuum on n^, having much smaller Ara^ or 
large oscillatory frequencies, will be observationally distinguishable from that caused by the 
phase transition in IR DBI inflation. In addition, there are no observable non-Gaussianities 
associated with them. 

We next look at the effect of the non-BD vacuum on the IR DBI inflation. Now e is 
negligibly small and s = The natural scale of M is the red-shifted string scale, so 

M oc 1/iVf^^ and = As we see from Sec. [3A1 iV°^^ is typically smaller than Ne. 

So the variation of Us can be comparable to (9(0.1). Using ( 1C.14I) . over (9(10) e-folds the 
change of the arguments in the trigonometric function is 2Mcs/H ^ 1, so the modulation is 
oscillating rapidly. In fact since the larger scales are associated with smaller M during infla- 
tion, by the time that the modulation stops oscillating, we have M ~ H/cs- This is already 
beyond the validity region of the non-BD vacuum calculations, and in fact is the place where 
the stringy phase transition takes effect. So we conclude that, in IR DBI model, the effect 
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of the non-BD vacuum is smoothly connected to the phase transition. Phenomenologically 
they look different from the phase transition by causing frequent oscillations in the spectral 
index. In Fig. [13] we illustrate this effect on the power spectrum. The amplitude of the 
oscillatory modulation increases and the frequency decreases towards the large scales, and 
finally merges into the phase transition. 

The effect of the non-BD vacuum on the large non-Gaussianities in DBI inflation is 
studied in [65]. The distinctive signature is the rising behavior in the shape of the 3- 
point function in the folded triangle limit. This may also be an interesting clue to a better 
understanding of the properties of the non-Gaussianity during the phase transition. 



D Details of numerical calculations 

The zero mode motion of the brane in a warped throat is captured by the DBl-CS action 
(13.11) . Varying the action, the exact form of the equations of motion is given by, 

3/' 2 /i7;v 3 \ , , /' 1 / , f 



V' + ^ , (D.l) 



Hn = -\lHct>l, (D.2) 

f{^) = ^, X^nnnR^. (D.4) 

In the above differential equations, we choose to use the number of e-folds = Ina(t) as the 
time coordinate, as the intrinsic time scale of inflationary dynamics is the Hubble time H~^. 
The subscript at denotes derivatives with respect to A^e, i-^- (t'N = dcp/dNe, (pNN = d'^(p/dN^. 
We also denote derivatives with respect to </> by a prime, i.e. f = df/d(f), V = dV/dcf). 

The equations flDlj) and can be integrated numerically using the conventional 



Runge-Kutta method. We note that ignoring Eq. flD.2p and setting if^r /i? = in Eq. fID.ip 
does not introduce detectable errors to the results of numerical calculation, because H can 
safely be treated as a constant for IR DBI inflation. Nevertheless, we have put Eq. (1D.2P 
through numerical integration together with (ID.ip for self-consistency. 



To integrate the equation of motion, the code needs to know the initial values 0(0) and 
07V (0). Due to the attractor behavior of the IR DBI dynamics, these initial conditions 
will be irrelevant to our calculation of observables as long as inflation lasts a few e-folds 
more than the minimum number required to solve the horizon problem. This can always 
be done in the IR DBI model, since we have the freedom to extend the start of inflation 
to the IR end of the throat by choosing 0(0) appropriately. In practice, we choose 0(0) < 
i^A/Ae/SO, so that according to (13. 3p . we will have at least 80 DBI e-folds, and roughly 90 
total e-folds (assuming 10 non-relativistic e-folds). This is good enough to make sure that, 
when we calculate primordial density perturbations on scales relevant to CMB temperature 
anisotropies (roughly 50 ~ 60 e-folds before the end of inflation), the inflationary background 
is well on the attractor solution. 
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To set up the model, the numerical code needs the five microscopic parameters {Nb, 
ub, nAh\, Qs ^^'^rus, P}. The five input parameters need to satisfy various bounds for model 
building consistency. The following microphysical bounds are imposed in the code: 

• the geometric constraint from compactification (Eq. (13. 32^ ) 




• the maximum number of branes generated by antibrane-fiux annihilation (Eq. (13.311) ) 

• the lower bound on string scale 

msgj^^^ > TeV ; 

• the lower bound on infiation scale 

nAh\T, > TeV^ ; 

• the warp factor Ha < 1, 

nAh\ < Ub . 

(In these bounds, we fix the string coupling Qs = 0.1 and = 1. The effects of different Qs 
and ttB are discussed in Sec. 13.51 ) Before calculating the density perturbation, the code per- 
forms checks on all of the above bounds to make sure the input parameters are theoretically 
consistent. 

After numerically integrating the inflaton equation of motion, our immediate result is 
0(iVe), 7(-/Ve) and H{Ne). Then we use (I3.28P to calculate the curvature perturbation Pfi{Ne) 
generated during inflation. Once we have PR(iVe), the horizon crossing relation (13.301) can 
translate P^^Ne) to PR^k), (noting Ne = Ntot—Ne,) which is then fed into CAMB to generate 
the CMB temperature anisotropy spectrum and the matter power spectrum. 
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